The coexistence of superconductivity and ferromagnetism in nano-scale metallic grains 
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A nano-scale metallic grain in which the single-particle dynamics are chaotic is described by the 
so-called universal Hamiltonian. This Hamiltonian includes a superconducting pairing term and 
a ferromagnetic exchange term that compete with each other: pairing correlations favor minimal 
ground-state spin, while the exchange interaction favors maximal spin polarization. Of particu- 
lar interest is the fluctuation-dominated regime where the bulk pairing gap is comparable to or 
smaller than the single-particle mean level spacing and the Bardeen-Cooper-Schrieffer theory of 
superconductivity breaks down. Superconductivity and ferromagnetism can coexist in this regime. 
We identify signatures of the competition between superconductivity and ferromagnetism in a num- 
ber of quantities: ground-state spin, conductance fluctuations when the grain is weakly coupled 
to external leads and the thermodynamic properties of the grain, such as heat capacity and spin 
susceptibility. 

PACS numbers: 74.20.-z, 75.75.-c, 74.78. Na, 74.25. Bt 
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I. INTRODUCTION 

Advances in nanofabrication techniques in the mid- 
1990s have made it possible to study the properties of in- 
dividual ultra-small metallic grains with linear dimension 
of several nanometers (see Ref. [l| for a review). Discrete 
energy levels of such nano-scale grains were measured 
by single-electron tunneling spectroscopy [HI]. Grains 
made of various materials have been studied, probing dif- 
ferent regimes of electron-electron interaction and spin- 
orbit scattering. Recent experimental techniques allow 
for better control of the shape and size of the grain [j| . 

Superconductivity in bulk metals was explained by the 
Bardeen-Cooper-Schrieffer (BCS) mean- field theory Q. 
Superconducting correlations in a finite-size system are 
characterized by the ratio of the bulk pairing gap A to 
the mean single-particle level spacing S. BCS theory is 
valid in the bulk limit A/S ^> 1 but breaks down in 
the fluctuation-dominated regime A/S < 1. In the first 
spectroscopy experiments on metallic grains @-Bj], both 
regimes were accessed by studying aluminum grains of 
different sizes. A pairing gap was observed in the excita- 
tion spectra of the larger grains with an even number of 
electrons. However, in the smallest grains, such a pair- 
ing gap could not be resolved. These experiments gen- 
erated much interest in superconductivity in finite-size 
systems and, in particular, in the fluctuation-dominated 
regime. It was proposed that pairing correlations in the 
fluctuation-dominated regime can be observed through 
the number-parity dependence of thermodynamic observ- 
ables such as the heat capacity and spin susceptibility @- 
[To| . For example, the spin susceptibility of an odd grain 
exhibits a re-entrant behavior (i.e., a local minimum) as 
the temperature decreases 0). 

The linear size of the ultra-small grains studied in spec- 
troscopy experiments is typically smaller than the mean 
free path of the electrons inside the grain, and the single- 
particle dynamics are dominated by scattering from the 



boundaries of the grain. In grains with irregular bound- 
aries, the single-particle dynamics are chaotic. A charac- 
teristic energy scale in such grains is the Thouless energy 
i?Th, which is determined by the time required for an 
electron at the Fermi energy to cross the grain. Energy 
spectra of chaotic grains fluctuate from sample to sample. 
In an energy window of width i?Th around the Fermi en- 
ergy, the statistical properties of the single-particle ener- 
gies and wave functions are described by random-matrix 
theory (RMT) [TMl- 

When the dimensionless Thouless conductance gxh = 
ETh/S is large, a metallic grain with chaotic single- 
particle dynamics is described by the so-called univer- 
sal Hamiltonian [l5|, [l6[ . The single-particle part of this 
Hamiltonian follows RMT, while the interaction is com- 
posed of three universal terms: a classical charging en- 
ergy term that depends on the capacitance of the grain, 
a BCS-like pairing term (which is suppressed in the pres- 
ence of an external orbital magnetic field), and a spin 
exchange term (in the absence of spin-orbit scattering). 

The single-particle part of the universal Hamiltonian 
and the pairing term favor minimization of the total spin, 
while the ferromagnetic exchange term tends to polarize 
the grain. In the bulk limit A/S 3> 1, the ground state 
for an even particle number is either completely paired or 
completely polarized, depending on the ratio J s /S, where 
J s is the spin exchange coupling constant. There is a 
first-order phase transition between the superconducting 
and ferromagnetic phases, whose signature is the macro- 
scopic onset of magnetization due to the Stoner insta- 
bility. In a finite system, there is strictly speaking no 
phase transition, and the ground-state spin of the grain 
increases stepwise as a function of J s /S In partic- 

ular, in the fluctuation-dominated regime A/S < 1, there 
is a region in the ground-state phase diagram where the 
electrons in the grain are partly paired and partly polar- 
ized, i.e., superconducting and ferromagnetic correlations 
coexist in the ground-state wave function (T3, HH ■ 
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Here we discuss various signatures of the competition 
between pairing and exchange interactions in ultra-small 
metallic grains that are described by the universal Hamil- 
tonian in the absence of orbital magnetic field and spin- 
orbit scattering fl?l fT^42lj. 

The outline of this paper is as follows. In Sec. [TT] we 
discuss the universal Hamiltonian and its general solu- 
tion. In Sec. Mil we describe the ground-state phase dia- 
gram of the universal Hamiltonian with an equally spaced 
single-particle spectrum [TtJ and discuss mesoscopic fluc- 
tuations of the ground-state spin. In Sec. |lV]we discuss 
signatures of the competition between superconductivity 
and ferromagnetism in the mesoscopic fluctuations of the 
linear conductance of a grain that is weakly coupled to 
external leads [l9j. In Sec. |Vj we discuss the signatures 
of this competition in thermodynamic observables of a 
nano-size grain and the mesoscopic fluctuations of these 
observables HIHD]. We conclude in Sec. EH 



to the interaction terms of the universal Hamiltonian (fT]) 
(where we have omitted the charging energy term re- 
lated to vq). The fluctuating part of the interaction is 
suppressed by 1 / grh [HI, EH and can be ignored in the 
limit <7Th > !• 

In practical calculations, the computational effort can 
be reduced by truncating the band width of the single- 
particle space. The pairin g co upling constant G is then 
renormalized according to [23l . l24j 
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where iV S p is the number of single-particle orbitals. The 
exchange coupling constant J s /S (expressed in units of 
the mean level spacing), is a material-dependent param- 
eter. Its values for different elements were tabulated in 
Ref. H. 



II. THE UNIVERSAL HAMILTONIAN 

In the absence of spin-orbit scattering and orbital mag- 
netic field, the universal Hamiltonian for a fixed number 
of electrons has the form [III, EH 

H= 12 (1) 

i,<7=t,4. 

where 

^ t= E c lt c 4 and ^-E c ^t ( 2 ) 

i i 

are the pair creation and annihilation operators, and S 
is the total spin of the electrons in the grain. The single- 
particle energies Cj follow the statistics of the Gaussian 
orthogonal ensemble (GOE) of RMT p]]. When the par- 
ticle number N is not fixed (e.g., in the presence of a 
varying gate voltage), an additional charging energy term 
e 2 N 2 /2C (where C is the capacitance of the grain) must 
be included in (UJ). 

The randomness of the single-particle eigenstates \i) 
induce fluctuations in the electron-electron interaction 
matrix elements Vhm, and the latter form an induced 
two-body ensemble (22J. The universal terms in (fT]) fol- 
low from the average part of this induced random inter- 
action. This can be derived from the requirement that 
the average interaction Vij-ki be invariant under a change 
of the single-particle basis (the GOE, which describes the 
one-body part of the Hamiltonian, is invariant under an 
orthogonal transformation of the single-particle basis). 
There are three orthogonal invariants: S^Sji, SijSki and 
SuSjk- We then have 

%;M = voSikdji - GSijSki + JsSuSjk , (3) 

where Vq,G and J s are constants. When expressed in 
its second-quantized form, this average interaction leads 



A. Solution to the universal Hamiltonian 

The universal Hamiltonian is an integrable model, 
which can be solved by generali zing Richardson's solu- 
tion of the pairing proble m 12a . |27| to include the ex- 
change interaction (T3, HH, l29l|. Below we discuss this 
solution and the quantum numbers used to classify the 
many-body eigenstates. 

We decompose the set of single-particle orbitals into a 
set B of singly occupied orbitals and the complementary 
set U of empty and doubly occupied orbitals. The pairing 
interaction scatters only time-reversed pairs of electrons 
from doubly occupied orbitals to empty orbitals. We can 
thus solve a reduced BCS-like Hamiltonian within the set 
U 

Hu= eiC l c i* ~ Gp ^P ■ ( 5 ) 

We denote the eignestates of |5| by \Q)u- The singly 
occupied orbitals in B are not affected by pair scattering 
and are referred to as "blocked" orbitals. They thus form 
a set of good quantum numbers, and we can solve for a 
reduced spin exchange Hamiltonian within the set B 

Hb = 22 ticLcia ~ JsS 2 ■ (6) 

The eigenstates of © are found by coupling the spins 
1/2 of the electrons in the set B to total spin S and spin 
projection S z . We denote the corresponding eigenstates 
by \jSS z )b, where 7 are quantum number used to dis- 
tinguish between states with the same SandS 2 [H[l|. 
The orbitals in U have spin zero and do not contribute 
to the total spin. The Hamiltonians in ([5]) and (|6]) com- 
mute, and the eigenstates of the universal Hamiltonian 
(TTJ are a direct product \()u ® \"fSS z )jg with energies 
E = Eu + Eg . 



The eigenstates of the reduced BCS model ([5]) can be 
found by using Richardson's solution [26|, [13] ■ A solu- 
tion with m pairs is expressed in terms of m Richardson 
parameters i? M that satisfy the following m coupled non- 
linear equations (jj, = 1, . . . m) 

G + 2 51 R u -R u = ^ 2e, - i? M ' ^ 

The eigenenergy of the reduced BCS-like Hamiltonian ([5]) 
is the sum of these parameters 

m 

Eu = ^ Rfi , (8) 
and the corresponding eigenstate can be written as 

io« « n (e 2^ c ^ c |o) • (9) 

In the limit G — > 0, different configurations of doubly 
occupied orbitals are not mixed, and for any Richardson 
parameter there is a doubly occupied orbital k such 
that — > 2efc. There is a one-to-one correspondence 
between this non-interacting solution and a solution for 
a finite pairing coupling constant G. Thus the quantum 
numbers £ of a pairing solution can be chosen to be the 
set of doubly occupied orbitals in the corresponding non- 
interacting solution. In fact, the values of R^ can be 
found numerically by evolving the initial non-interacting 
values as a function of increasing G. 

For a given set B with b — N — 2m blocked orbitals, 
the number of different eigenstates \jSS z )is is 2 b . Their 
energies 

E B = J2^-JsS(S + l) (10) 

are independent of_7. The number of such states with 
given 5 and S z is [28j 

db ^=(s + b/2)-(s + l + b/2) ■ ( U ) 

III. GROUND-STATE SPIN 

In this section we present results for the ground-state 
spin of the universal Hamiltonian (JIJ as a function of 
the pairing gap A/5 and exchange interaction strength 
J s /S (both measured in units of 5). We find the lowest 
energy E(S) for a given spin 5 using the method dis- 
cussed in Sec. Ill Al and then minimize E(S) with respect 
to 5. The ground-state spin of the grain is determined 
by the competition between various terms in the univer- 
sal Hamiltonian. The one-body part (kinetic energy plus 
one-body confining potential) and the pairing interaction 
favor minimal spin, while the exchange interaction favors 
a maximally polarized state. 
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FIG. 1: Ground-state phase diagram of a grain with an 
equally spaced single-particle spectrum in the Js/S- A/ 5 plane 
for an even (left panel) and odd (right panel) number of elec- 
trons. The lines separate sectors of different spin values (the 
numbers are the corresponding spin values). We observe an 
intermediate regime (SC-FM) in which the ground state of 
the grain is partly polarized and partly paired. Reproduced 
from Ref. 



A. Equidistant single-particle spectrum 

Results for a generic equidistant single-particle spec- 
trum with level spacing 8 are shown in Fig. [TJ We find 
three different phases: a superconducting phase where 
the number of pairs is maximal, a ferromagnetic phase 
where the system is fully polarized with S = N/2, and 
an intermediate regime where exchange and pairing cor- 
relations coexist. The coexistence regime describes a 
partially polarized state, in which 2S electrons reside in 
singly occupied levels closest to the Fermi energy and the 
remaining electrons are paired to give spin zero [13, EH ■ 
The existence of pairing correlations in this intermediate 
regime is reflected in the shift of the spin transition lines 
to higher values of the exchange interaction strength as 
the pairing gap A/5 is increased. Fig. Q] demonstrates the 
stepwise increase of the ground-state spin from its min- 
imal value S = (even grain) or 1/2 (odd grain) to its 
maximal value of S = N/2 as we increase the exchange 
coupling constant J s /S at fixed A/5. 

In the absence of pairing correlations (i.e., A = 0) and 
for an equidistant single-particle spectrum, the transi- 
tion from spin S to spin 5+1 occurs for an exchange 
coupling of J s /S = (25 + l)/(25 + 2). This stepwise in- 
crease is known as the mesoscopic Stoner staircase [l5[ 
and is shown in Fig. [2] An interesting qualitative change 
in the presence of pairing correlations is the possibility 
of a spin jump (i.e., A5 > 1) for the first step. All sub- 
sequent steps have A5 = 1. For A/5 < 0.6, the first step 
corresponds to A5 = 1. However, for 0.6 < A/5 < 0.8, 
the ground-state spin changes from 5 = to 5 = 2 in 
one step for a value of the exchange coupling in the range 
0.87 < J s /S < 0.9. The height of the first step gets larger 
with increasing A/5. For example, in Fig. [2] the spin in- 
creases from 5 = to 5 = 3 in the first step at A/5 = 1 
and then continues to increase by one unit at a time. 




FIG. 2: Ground-state spin versus the exchange coupling J a /8 
for an even grain with an equally spaced single-particle spec- 
trum at A/S = 0, 1,2. As A/S increases, the staircase func- 
tion gets more compressed and pushed to higher values of 
Js/5. It also exhibits spin jumps, e.g., AS = 3 for A/8 — 1. 
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FIG. 3: Ground-state spin distribution Ps for a metallic grain 
with different values of A/ 5 and Js/S. The darker (lighter) 
histograms describe the even (odd) grains. 



B. Mesoscopic fluctuations 

The randomness of the single-particle part of the uni- 
versal Hamiltonian leads to mesoscopic fluctuations of 
the ground-state spin. Features of the resulting ground- 
state spin distribution Ps were studied in Ref. [30]. 

Fig. [3] shows the ground-state spin distribution for sev- 
eral values of the parameters A/S and J s /S. In the ab- 
sence of interactions, the ground-state spin will always 
be minimal because of the cost to promote an electron 
out of the Fermi level to the next unoccupied level (to 
form a larger spin state). This energy cost undergoes 
mesoscopic fluctuations (because of the GOE nature of 
the single-particle spectrum), and at finite exchange in- 
teraction can sometimes be compensated by the gain in 
exchange energy (as the spin increases). Grains with 
single-particle levels close to the Fermi energy can ac- 
quire a non-minimal ground-state spin. However, in 
grains with no levels sufficiently close to the Fermi en- 
ergy, the gain in exchange energy cannot compensate 
the cost to promote electrons to higher levels, and the 
ground-state spin remains minimal. As a result the spin 
distribution Ps acquires a nonzero variance. This vari- 
ance increases with the exchange coupling constant be- 
cause higher ground-state spin values become more prob- 
able. The pairing interaction counteracts this behavior 
and favors smaller values for the ground-state spin, i.e., it 
suppresses the probability for large spin values. This can 
be clearly seen when comparing the left (A/ 5 = 0) and 
right (A/5 = 0.25) columns in Fig. [3] Similar results for 
the mesoscopic ground-state spin distribution were found 
in Ref. 0. 

A particularly interesting signature of the mesoscopic 
interplay of ferromagnetism and superconductivity in 
chaotic metallic grains is shown in Fig.|4j In this figure we 
show the spin distribution at a relatively large exchange 
interaction strength J s /5 — 0.9 (close to the Stoner in- 
stability) and a moderate pairing strength A/S = 1. This 



case is within the coexistence regime in the ground-state 
phase diagram for an equidistant spectrum (see Fig. [T]) 
and in the vicinity of a large spin jump. We observe in 
Fig. Q] that the even-grain spin distribution has a local 
maximum at S = 3 (besides the maximum at S — 0) and 
a dip at S — 1. Approximately 43% of the grains have 
spin S = and 5% have spin S — 1, but 9% have spin 
S = 3. Thus the probability for certain spin values is 
anomalously suppressed. We interpret such a dip in the 
ground-state spin distribution as a mesocopic signature 
of spin jumps. Such a dip has a width of about AS" ~ l 
and occurs only at values of the exchange interaction 
strength that are close to the Stoner instability. 

IV. TRANSPORT THROUGH A WEAKLY 
COUPLED GRAIN 

In this section we investigate the signatures of the 
coexistence of superconductivity and ferromagnetism in 
experimentally measurable quantities, i.e., the tunnel- 
ing conductance of a grain that is coupled to external 
leads. We assume weak coupling between the grain and 
the leads and consider the sequential tunneling regime, 
where 8 and the thermal energy associated with temper- 
ature T are much larger than a typical tunneling width 
r, i.e, 8,T HT (here and in the following we set the 
Boltzmann constant ks — 1)- In this regime the coher- 
ence between the leads and the grain can be ignored, and 
the conductance can be calculated using a rate equation 
approach (3ll . [32j. The input in this approach are the 
many-particle energies and wave functions of the univer- 
sal Hamiltonian, and here we obtain them using an ex- 
act diagonalization method. When the charging energy 
is much larger than the thermal energy, i.e., e 2 /C T, 
the conductance displays a series of sharp peaks as a 
function of gate voltage. This is the so-called Coulomb 
blockade regime (for reviews see, e.g., Refs. 0). We 
studied the statistics of these conductance peaks in the 
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FIG. 4: Ground-state spin distribution Ps for a metallic grain 
with A/8 — 1 and J„/8 — 0.9. Results are shown for both 
even (darker histograms) and odd (lighter histograms) parti- 
cle number. Spin jumps observed in the ground-state phase 
diagram of Fig. [1] (for a grain with an equidistant single- 
particle spectrum) manifest as a dip in the mesoscopic prob- 
ability distribution. 



quantum regime at a typical experimental temperature 
of T ~ 0.1 8. For each realization of the one-body Hamil- 
tonian, we calculated the five lowest many-particle eigen- 
states of the Hamiltonian ([1} using the Lanczos method 
for a bandwidth of iV S p = 17 single-particle levels and up 
to N — 19 electrons. For each tunneling event, the max- 
imal value of the linear conductance as a function of gate 
voltage is found numerically. We used ~ 4000 samples to 
collect sufficiently good statistics. 



A. Conductance peak spacings 

When the total energy of a grain is approximated 
by the classical charging energy e 2 N 2 /2C, the Coulomb 
blockade peaks are equally spaced with a spacing of e 2 /C. 
The actual spacing fluctuates and is given by e 2 /C + A2, 
where A 2 incorporates the effects of the discrete single- 
particle spectrum and interactions beyond the classical 
charging energy. The quantity A2 is known as the peak 
spacing and below we discuss its statistics. 

The distributions of A2 are shown in Fig. [5] for various 
values of A/5 and J s /5. For A = J s = and low temper- 
atures T <C 8, the peak spacing distribution is bimodal 
because of a number-parity effect caused by the spin de- 
generacy of the single-particle levels. The exchange in- 
teraction suppresses this bimodality [33, 34] since it leads 
to fluctuating spin polarization as discussed in Sec. lIIIB] 
At J s /8 — 0.6, the bimodality is completely washed out. 

The most significant effect of pairing correlations on 
the peak-spacing distribution is to restore this bimodal- 
ity. For an exchange coupling of J s /8 = 0.6, bimodal- 
ity starts to reappear for A/8 = 0.5. At even stronger 
pairing A/ 8 = 1, the peak spacing distribution is well 
separated into two components. The left part describes 
the sequence of even-odd-even (EOE) tunneling events, 
while the right part corresponds to odd-even-odd (OEO) 
transitions. For A/8 > 1 the ground-state spin of an 
even (odd) grain is almost always S = (S = 1/2), and 




FIG. 5: Peak-spacing distributions obtained from exact di- 
agonalization at T/8 = 0.1 for several values of A/5 shown 
in the panels. For each value of A/S we show results for 
J s = (dashed, grey-filled histograms) and J s /5 = 0.6 (solid 
histograms). For A = we also compare with the analytic 
result at T <C <5 and J s = (solid line) [14| . Reproduced from 
Ref. fijl. 



mesoscopic fluctuations of the spin can be ignored. In- 
creasing the exchange coupling constant merely shifts the 
two components of the distribution closer together. The 
separation of the peak-spacing distribution into two parts 
in the presence of pairing correlations can be understood 
qualitatively using a fixed-S* BCS mean-field approxima- 
tion [l|. In an EOE sequence, the first peak in the con- 
ductance corresponds to the blocking of an additional 
single-particle level and the second peak to the removal 
of this blocked level by creating an additional Cooper 
pair. In an OEO sequence, these two events are reversed, 
leading to a separation of the peak spacings by about 
8A 2 » 4A - 3J S (assuming A > 8, J s ) pj|. The contri- 
bution of the exchange energy is 3 J s since in this regime 
of strong pairing correlations the ground-state spin is al- 
most always S — (S = 1/2) for the even (odd) grain. 



B. Conductance peak heights 

Distributions of conductance peak heights g max are 
shown in Fig. [5] For A = J s = and very low tem- 
peratures T -C 8, the peak height distribution follows 
from the Porter-Thomas distribution of eigenvector com- 
ponents of a random matrix and is known analytically 
[35j. The exchange interaction increases the number of 
states contributing to the conductance by bringing down 
higher spin states. This results in a narrower peak height 
distribution. 

The mean value g max and the standard deviation (rel- 
ative to the mean) cr(<? max )/<7 max of the peak height g max 
are shown in Fig. [7] For small pairing strengths both 
quantities are strongly suppressed by the exchange in- 




ln(g /g ) 

FIG. 6: Peak-height distributions obtained from exact di- 
agonalization at T = 0.15 for J s — (dashed, grey- filled 
histograms) and J s /S — 0.6 (solid histograms). Different 
panels describe different values of A/5. We compare with 
the analytic result at T <C 8 and J s = (solid line) (35[. 
Since <? max fluctuates over several orders of magnitude, we use 
ln(g ma x/gmax) where (? max is the average conductance peak 
height. Reproduced from Ref. [lSj. 




FIG. 7: Average (left panel) and relative width (right 
panel) [l^| of the peak height distribution as a function of 
A/5 at T/5 = 0.1 for J s = (open circles) and J s /8 = 0.6 
(solid circles) obtained from exact diagonalization. In the left 
panel we compare with the results of an equidistant single- 
particle spectrum (solid lines). The vertical bars in the right 
panel show the sampling error of the GOE distribution with 
4000 samples (the errors in the left panel are smaller than the 
symbols). The right panel is adapted from 



is known to be superconducting in granular form 37 1 
and has a relatively large exchange interaction strength 
of J s /5 ~ 0.59 - 0.72 0. 



teraction. This suppression was observed in semicon- 
ductor quantum dots, where the Cooper pairing channel 
can be ignored and a closed solution for the conductance 
is available [28|. There, the finite-temperature suppres- 
sion of the probability of small conductance peak heights 
induced by exchange correlations led to a good agree- 
ment between theory and experiment at low tempera- 
tures m m . 

The effect of pairing correlations on the excitation 
spectrum of the grain is to induce a gap that pushes states 
with large spins to higher energies. Thus, we observe that 
for J s /6 = 0.6 the probability of small conductance peak 
heights and the standard deviation of the peak height 
distribution increase with A/ 5 up to A/ 5 ~ 1. In par- 
ticular, if the pairing interaction is sufficiently strong to 
destroy all spin polarization, i.e., at A/5 = 1.0, the peak 
height distribution becomes independent of the exchange 
interaction strength J s /5 (see Fig. [B]). 

We define a mesoscopic coexistence regime of super- 
conductivity and ferromagnetism by the simultaneous oc- 
currence of a bimodality in the peak spacing distribution 
(a signature of pairing correlations) and the suppression 
of peak height fluctuations (a signature of ferromagnetic 
correlations). In Figs. [51 [5] and [7] we observe such a meso- 
scopic coexistence case for A/6 = 0.5 and J s /5 = 0.6 at 
a typical experimental temperature of T = 0.1 6. Inter- 
estingly, for an equidistant spectrum this grain still be- 
longs to the superconducting regime (see Fig. [1]). We 
conclude that mesoscopic fluctuations enlarge the pa- 
rameter regime for which coexistence of pairing and ex- 
change correlations can be measured. A good material 
for observing the coexistence regime is platinum, which 



V. THERMODYNAMICS 

Here we discuss signatures of the interplay between su- 
perconductivity and ferromagnetism in thermodynamic 
observables of the grain. 

Effects of pairing correlations on the heat capacity 
and spin susceptibility of ultra-small metallic grains have 
been studied theoretically using a number of methods @- 
ITol . 1241 [38M40I ] . Thermodynamic observables of small 
metallic clusters were studied experimentally in Ref. [4l[ . 
Signatures of pairing correlations in the bulk limit A/S ^> 
1 are well described by the grand-canonical BCS theory. 
These signatures are the exponential suppression of both 
the heat capacity and spin susceptibility at low temper- 
atures because of the gap in the density of states, and 
a sharp peak in the heat capacity near the critical tem- 
perature T c , which is a signature of a second-order phase 
transition (see the BCS results in Fig. [5]). 

In a small isolated grain, the fluctuations of the pair- 
ing gap and the dependence of observables on the parity 
of particle number (he., odd-even effects) become impor- 
tant 0-ES Hi H, EH- An example is shown in Fig. |H 
where the heat capacity (top row) and spin susceptibil- 
ity (bottom row) for the pairing model with an equally 
spaced single-particle spectrum are shown versus temper- 
ature for A/6 = 0.5 (left column) and A/5 = 3.0 (right 
column) . The results for even (odd) grains are shown by 
solid (dashed) lines, while the grand-canonical BCS re- 
sults are shown by the dash-dotted lines (for A/5 = 3.0). 
Finite-size effects result in smoothing of the sharp peak 
in the even-grain heat capacity and its gradual transfor- 
mation into a shoulder as A/5 decreases. This shoulder 
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FIG. 8: Heat capacity (top row) and spin susceptibility (bot- 
tom row) versus temperature. The spin susceptibility is mea- 
sured in units of the Pauli susceptibility \p — 2fi%/8, where 
(iB is the Bohr magneton. Shown are results for supercon- 
ducting grains with an equally spaced single-particle spec- 
trum, and with A/8 — 0.5 (left column) and A/8 = 3.0 (right 
column). There are no exchange correlations (i.e., J s = 0). 
Solid and dashed lines describe grains with even and odd par- 
ticle number, respectively. The observables are calculated us- 
ing Richardson's solution [H|,[27j at low temperatures and the 
method of Ref. [2fl ] at higher temperatures. The dash-dotted 
lines (for A/ 8 = 3) are the grand-canonical BCS results. 



a pure pairing Hamiltonian (i.e., in the absence of an 
exchange interaction). 

Below we discuss this semi-analytic method and 
present the main results for the heat capacity and spin 
susceptibility of the grain [2l| . 

A. Spin and number-parity projections 

The grand-canonical trace of an operator can be writ- 
ten as a sum of diagonal matrix elements of the operator 

over a complete set of quantum numbers Ai, A2, We 

define Tr.\ to be a trace that is restricted to the partic- 
ular values of a subset A of this complete set (while all 
other quantum numbers are summed over). An impor- 
tant example is when A is the spin 5* of the grain. For 
an operator X we can write TrA = ^2 S TrgA when Trs 
denotes a trace over all states with fixed total spin S. 
The following identity allows us to express the trace of a 
scalar operator X at fixed spin 5* in terms of the traces 
at fixed values of the spin component S z 

Tr s A = (2S + I)(Trs z= sA-Trs 2=s+1 A) . (12) 

A similar identity applies to canonical traces at fixed par- 
ticle number N and spin S. 

Using the identity (TT2"j) , we can write the canonical 
partition function Zn(J s ) of the universal Hamiltonian 
(Q} in the form 



disappears when A/ 5 < 1, but the heat capacity con- 
tinues to be enhanced for an even grain (as compared 
with an odd grain) at intermediate temperatures. The 
unpaired electron in the odd grain leads to a Curie-like 
divergence of the spin susceptibility at low temperatures 
(x ~ l/^ 1 )- The combination of this divergence and the 
suppression of x due to pairing correlations produces a 
re-entrant behavior (i.e., a local minimum) of the odd- 
grain spin susceptibility as a function of temperature Q ■ 
This signature of pairing correlations exists even in the 
fluctuation-dominated regime A/S < I. 

Exchange correlations affect the thermodynamic ob- 
servables shown in Fig. [H Signatures of the competi- 
tion between exchange and pairing correlations in ther- 
modynamic observables were studied in Refs. [2(| and 
(2lj . The case of an equally spaced single-particle spec- 
trum was studied in [20| using a quantum Monte-Carlo 
method, while the general case of a fluctuating RMT 
single-particle spectrum was considered in plj using a 
semi- analytical method. The latter method is accurate 
and efficient, making it useful for studying the mesoscopic 
fluctuations of thermodynamic observables. It uses a 
spin projection technique [28|, |42j to treat exactly the 
exchange interaction in ([T]) , and a functional integral ap- 
proach to treat approximately the pairing part of ((T|). 
The odd-even effect is reproduced by a number- parity 
projection. A similar approach was used in Ref. [lOj to 
study odd-even effects in thermodynamic observables for 



Z N (J S ) =^(25 + 1) 
s 

x e^ s < s+1 ) (Z N , S . =S ~ Z NA =s+i) , (13) 
where (3 — 1 /T and 

Z N , Sz = Tr/v,s, (e-^ BCS ) (14) 

is the S* z -projected canonical partition function of the 
BCS-likc Hamiltonian 

H BC s = J2 ei4* c i<r-GPiP. (15) 

»,<7=f4 

A similar expression can be written for the spin suscep- 
tibility (ps is the Bohr magneton) 

X ^ ( ^^_A_£ S(S+1)(2S+1) 

x e^ s ( s+1 ) (Z„ t s B=s ~ Z N ^=s+i) ■ (16) 

The heat capacity can be calculated from the partition 
function by taking numerical derivatives with respect to 
temperature. Thus the calculation of the heat capacity 
and spin susceptibility of a metallic grain described by ([T]) 
reduces to the calculation of the ^-projected partition 
function (fT4"|) . 
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The trace at fixed S z can be calculated using a discrete 
Fourier transform for the S z projection operator 



1 



Smax 



2S n 



1 ^ 

in— — o ma > 



(17) 



Here S max is the maximal possible value of the spin, and 
4> m — 27rm/(25 max 4- 1) are quadrature points. 

The canonical projection on a fixed particle number 
can be done using a similar Fourier transform, but it leads 
to cumbersome expressions in the framework of the path 
integral formalism described below. Important odd-even 
effects can be described by the simpler number-parity 
projection [43l - [45j 



P„ 



1 



1 



iirN 



(18) 



where rj = 1 (77 = — 1) corresponds to even (odd) number 
of particles. The canonical partition function for N par- 
ticles at fixed 77 and S z is then evaluated by the saddle- 
point approximation of the particle-number projection 
integral 



J N,r],S z 



f3dfi 



V,S Z 



Here 



Z vA =Tr n>Sr e -^Bcs-^) 



(19) 



(20) 



is the grand-canonical partition function at fixed values 
of 77 and S z - 



describes a non-interacting Hamiltonian in an external 
time-dependent pairing field A(r) (0 < r < (3) and T 
denotes time ordering. The saddle-point approximation 
to this integral (without the number-parity and S z pro- 
jections) gives the BCS theory. 

The auxiliary field A(r) can be separated into its static 
and time-dependent parts by expanding in a Fourier se- 
ries 



A(r) = A 



(23) 



where ui r = 2ttt / P (r integer) are the bosonic Matsubara 
frequencies. In the static-path approximation (SPA) (38l . 
EH El, we ignore the time dependance of A (i.e., we take 
A r = for all r ^ 0) and carry out the integral over the 
static field Aq exactly. In this approximation 



J v,S, 



/3rf|A c 
G 



^ /G)|Aol X5,(A ), (24) 



where Z^^JAo) = Tr Vi s z e~ l3HA ° is a static partition 
function. Using the explicit expressions (TIT)) and (|18p for 
the projection operators, and working in a Ao-dependent 
quasiparticle basis that diagonalizes H a , we find f2l"j 




B. Auxiliary-field path-integral formalism 

Here we present only the relevant definitions of the 
functional-integral formalism and summarize the results 

For more details see 



of the calculation of Z 



N,r),S z - 

Ref. [21|. 

We first discuss the evaluation of the grand-canonical 
partition function Z r/t s z in fl20p. Using the Hubbard- 
Stratonovich (HS) transformation 0, [47J , this partition 
function can be written as a functional integral over a 
complex auxiliary field A(r): 



J V,Sz 



Here 



X>[A, A^Tr^ 



s z 



Te 



-/dr(|A(r)| 2 /G + if A(T) ) 



(21) 



G 



e% - M - 



\A\piir 



f) 2 



|A | are the quasiparti- 



where Ei = 
cle energies. 

The SPA can be improved by carrying out the integra- 
tion over A r (with r ^ 0) in the saddle-point approxima- 
tion for each value of Ao. This leads to the random-phase 
approximation (RPA) correction to the SPA [3, |50| - [55| . 
This correction amounts to multiplying Z n ^ z (Ao) in (124)) 
by the factor 



c^>o) = n 



2Ei 



sinh 



(26) 



Here ±f2j are the eigenvalues of the 2N sp x 2A sp RPA 
matrix (N sp is the number of single-particle orbitals) 



ZErfij - f /i(7i7j + 1) 



''-fiiiaj + 1) - 2EiS i: 



where 



(27) 



and 



= 1 91nZ„, Sz (Ao) 
Ji P dEi 



(29) 



Next we discuss the approximate evaluation of the 
canonical partition function Zjq ^ s x in (|19|l - Substitut- 
ing the HS transformation pip for Z^ i s z in the integrand 



of (|19| and exchanging the order of integrations, we eval- 
uate the integral over [x in a saddle-point approximation 
for each fixed value of Ao- This saddle-point integration 
is applied to the grand-canonical static free energy 



F(A ) = \A \ 2 /G-T\nZ(A ) 



where 



Z(A Q ) = Tre- fiA *° = Jj4e _ ^ (e< -^ cosh ; 



(30) 



(31) 



The remaining ratio Z^ Sz (A )C^(A )/Z(A ) is 
treated as a pre-factor that does not enter in the saddle- 
point integration. We find 



Z 



N,v,S z 



f3d\A \ /2tt 



G \p 



d 2 F 



dfi 2 



-1/2 



x e-W°^ 2 e-^ N Z vA (A Q ) C^(A ) (32) 



with 

d 2 F 



dfi 2 



E 



/3E t Ui - n - f ) 2 + |A | 2 sinh(/3 J E i ) 



2£„ 3 cosh 2 



(33) 

The saddle-point equation for /u has the usual form of the 
particle-number equation 



TV 



OF 
dfi 



E 



H- 2 



tanh 



(34) 



but its solution fx = (i(Aq) depends on the value of the 
static field. 

We have shown that this number-parity projected 
SPA+RPA method is quite accurate by comparing it 
with exact results. The latter were calculated from the 
many-body eigenvalues of (JIJ using Richardson's solu- 
tion generalized to include the exchange interaction (see 
Sec. Ill A[) . In Fig. |H] we show the heat capacity (top 
panel) and spin susceptibility (bottom panel) for a grain 
with equally spaced single-particle spectrum and with 
A/5 = 1, J s /8 = 0.5. The exact results (lines) are 
compared with the SPA results (open symbols) and the 
SPA+RPA results (solid symbols) for both even and odd 
grains. These results demonstrate (i) the importance of 
the RPA correction, and (ii) the good accuracy of the 
SPA+RPA calculations. We note that Richardson's solu- 
tion is not practical at high temperatures since the num- 
ber of required eigenvalues increases exponentially. The 
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FIG. 9: Heat capacity C (top panel) and spin susceptibility \ 
in units of the Pauli susceptibility \p (bottom panel) versus 
temperature T/8 of a metallic grain described by the universal 
Hamiltonian (0 with equally spaced single-particle levels, and 
with A/5 = 1.0, Js/5 = 0.5. Results are shown for both even 
particle number (solid line, circles) and odd particle number 
(dashed line, squares). Results obtained using the generalized 
Richardson's solution of Sec. Ill Al (lines) are compared with 
the SPA results (open symbols) and the SPA+RPA results 
(solid symbols). 



deviation of the exact results from the SPA+RPA results 
at the higher temperatures in Fig. [3] is because the ex- 
act many-body eigenvalues were calculated only below a 
cutoff of - 30 S. 

The SPA+RPA method is not applicable below a cer- 
tain critical temperature, when the gaussian integration 
over one of the fluctuations A r diverges for at least one 
value of Ao, making the RPA correction unstable. It was 
recently proposed [56| that this problem can be over- 
come by treating non-perturbatively a low-energy collec- 
tive mode. At low temperatures, one can also use the ex- 
act solution (the number of eigenvalues required to calcu- 
late thermal observables at sufficiently low temperatures 
is not large). 



C. Heat capacity 

The main results for the heat capacity C as a function 
of temperature are shown in Fig. [TU] for several values 
of A/6 and J s /6. The symbols with vertical bars are, 
respectively, the average values C and standard devia- 
tions 5C calculated for an ensemble of 1000 realizations 
of single-particle GOE spectra in ([T|). The lines are the 
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FIG. 10: The heat capacity C versus temperature T/8 for an 
even grain (solid lines, circles) and for an odd grain (dashed 
lines, squares). Results are shown for grains with A/ 8 = 0.5 
(left column), A/8 = 1.0 (middle column) and A/S — 3.0 
(right column), and with J s /8 = (top row), J 3 /S = 0.5 
(middle row) and J 3 /S — 0.85 (bottom row). The symbols 
and vertical bars describe, respectively, the average value C 
and standard deviation 8C of the heat capacity for a fluctu- 
ating RMT single-particle spectrum in (flj. The lines corre- 
spond to an equally spaced single-particle spectrum, and the 
dash-dotted lines are the grand-canonical BCS results (where 
applicable) . 



results for an equally spaced single-particle spectrum (ob- 
tained from the exact solution at low temperatures and 
from the SPA+RPA at higher temperatures). The num- 
ber of single-particle levels N sp in the model space was 
odd and varied between 31 and 61 (depending on tem- 
perature) and the particle number was either N sp — 1 or 
Ngp (depending on its parity). The interplay between 
odd-even effects and mesoscopic fluctuations is further 
demonstrated in Fig. [TTJ The sum SC cvcn + <5C dd of the 
standard deviations for the even and odd grains is shown 
by dash-dotted lines, while the difference C CV on — Codd 
between the average values for the even and odd grains 
is shown by solid lines. 

Both figures \W\ and QT] show that the exchange inter- 
action shifts the odd-even effect to lower temperatures. 
In addition, when A < <5, the exchange suppresses the 
odd-even effect. However, when A > 8, the exchange 
suppresses only the right side of the even-grain shoul- 
der and therefore transforms this shoulder into a peak. 
Overall, the effect of exchange is less prominent when 
A > S; higher values of J s /5 are required to produce a 
visible change, and those changes occur mostly at higher 
temperatures. 
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FIG. 11: The sum SC even + SC dd of standard deviations 
(dash-dotted lines) and the difference Ccven — C a dd between 
average values (solid lines) of the heat capacities of an even 
grain and an odd grain are shown versus temperature. 



These results are consistent with the phase diagram 
discussed in Sec. IIII1 The dependence of a thermody- 
namic observable on J s /S arises through the dependence 
of the many-particle energies on J s /S, which is deter- 
mined by the total spin S of the states. A typical exci- 
tation energy of a state with non-ground-state spin in- 
creases with A/S and thus the change in the relative 
weight of such a state decreases with increasing A/ 6. At 
larger values of A/S, this relative change becomes ap- 
preciable only at higher temperatures, resulting in the 
suppression of the right side of the even-grain shoulder 
but in a much weaker dependence on J s /S of its left side. 

Another interesting question is whether the odd-even 
effect in the heat capacity can be resolved despite the 
mesoscopic fluctuations in experiments in which the num- 
ber parity of electrons in the grain is unknown. The 
condition for resolving the odd-even effect in such exper- 
iments is that Covon — Codd be larger than <5C OV on+<5Codd, 
namely, when the solid line in Fig. [TT] is above the 
dash-dotted line. We observe that in larger grains with 
A/S > 1, there is always a substantial temperature inter- 
val in which this condition is satisfied. In the crossover 
regime A/ 5 ~ 1, the odd-even effect starts to be washed 
out by mesoscopic fluctuations for a moderate strength 
of the exchange interaction. In the smaller grains with 
A/S = 0.5, the observability of the odd-even effect is 
marginal already in the absence of exchange interaction. 
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FIG. 12: The spin susceptibility \ in units of xj b = Xf/(1 ~~ 
Js/8) versus temperature T/S for even and odd grains for the 
same values of A/ 8 and J a /5 as in Fig. [TD] Symbols and 
lines follow the same convention as in Fig. [TO] but for the spin 
susceptibility. 



D. Spin susceptibility 

The results for the spin susceptibility x are shown in 
Figs. [TJ] and Q2] for the same ensemble of Hamiltonians 
used for the heat capacity. In both figures, \ is nor- 
malized by the J s -dependent high-temperature limit of 
its average value x.J s — Xp/(1 — Js/8) US HH, where 
Xp = 2fi^/S is the Pauli susceptibility. Our notation fol- 
lows the same convention as in Figs. [TO] and QT] except 
that the solid lines in Fig. IT51 correspond to the odd-even 
difference of spin susceptibilities (rather than the even- 
odd difference taken for the heat capacity). 

As with the heat capacity, the exchange interaction can 
shift odd-even effect in the spin susceptibility to lower 
temperatures. When A/S < 1, the exchange can com- 
pletely suppress the re-entrant behavior of the odd-grain 
average spin susceptibility. However, when A/S > 1, the 
re-entrant behavior is enhanced by the exchange interac- 
tion. This effect is similar to the suppression of the right 
side of the even-grain shoulder in the heat capacity, but 
it is more prominent for the spin susceptibility, whose 
high-temperature value is strongly affected by exchange 
correlations. 

In general, exchange correlations enhance the meso- 
scopic fluctuations of x- These fluctuations become par- 
ticularly strong (as J s /S increases) in the fluctuations- 
dominated regime, even when x is measured in the units 
of the already enhanced average value XJ B - These strong 
fluctuations of x m t ne smaller grains are explained by 
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FIG. 13: The sum 5xev<m+<5Xodd of standard deviations (dash- 
dotted lines) and the difference Xodd — Xeven between average 
values (solid lines) of the spin susceptibilities for an odd grain 
and for an even grain versus temperature. The susceptibility 
X is measured in units of \J S — Xp/(1 — Js/S). 

the large dispersion of the magnetization of the grain or 
alternatively by the strong dependence of x on excita- 
tion energies of states with different spin values. Fluc- 
tuations in such excitation energies result in very large 
fluctuations of x m the smallest grains. 

VI. CONCLUSIONS 

We have discussed various signatures of the interplay 
between superconducting and ferromagnetic correlations 
in nano-scale metallic grains. In particular, we have con- 
sidered chaotic grains with a large dimensionless Thou- 
less conductance, which are well described by the uni- 
versal Hamiltonian. The non-interacting part of this 
Hamiltonian is fluctuating and is modeled by random- 
matrix theory, while the coupling constants of its inter- 
action part do not fluctuate. For a fixed particle num- 
ber, this universal interaction consists of pairing and ex- 
change terms describing interactions in the Cooper and 
spin channels, respectively. The single- particle part of 
the Hamiltonian and the pairing term favor minimiza- 
tion of the total spin and lead to superconductivity in 
the bulk, while the spin exchange term tends to polarize 
the grain and leads to ferromagnetism when the exchange 
coupling constant exceeds the Stoner instability thresh- 
old. We have identified signatures of the competition 
between pairing and exchange correlations in the follow- 
ing observables: (i) the ground-state spin of the grain; 
(ii) the mesoscopic fluctuations of the conductance of a 
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grain that is weakly coupled to leads, and (iii) thermo- 
dynamic observables such as the heat capacity and spin 
susceptibility. We summarize these signatures as follows. 

(i) Ground-state spin. The ground-state phase diagram 
of the grain exhibits a coexistence regime in which the 
electrons in the grain are partly paired and partly polar- 
ized. The ground-state spin at fixed A/ S is a monotonic 
stepwise function of J s /5. The values of J s /S at which 
such steps occur increase with A/S. When A/5 is suffi- 
ciently large, the change AS of the ground-state spin in 
the first step can be larger than 1 (this is known as a spin 
jump). 

The mesoscopic fluctuations of the ground-state spin 
acquire a finite variance in the presence of exchange 
correlations, and they are suppressed by the compet- 
ing pairing correlations. When both the exchange and 
pairing interactions are sufficiently strong (A/ 5 ~ 1 and 
J s /S ~ 0.9), the ground-state spin distribution has a lo- 
cal minimum, which is a signature of spin jumps. 

(ii) Linear conductance. In the absence of pairing and 
exchange correlations, the distribution of the Coulomb- 
blockade conductance peak spacings is expected to be bi- 
modal (because the single-particle energy levels are spin 
degenerate). In semiconductor quantum dots (in which 
the pairing interaction is repulsive and its effects can be 
ignored), this bimodality is suppressed by the exchange 
interaction. However, the competing attractive pairing 
correlations in metallic grains can restore the bimodality 
of the conductance peak spacing distribution. 

The number of many-particle energy levels that par- 
ticipate in the low-temperature conductance increases in 
the presence of exchange correlations since levels with 
higher spin values shift down in energy. This leads to the 
suppression of the peak height fluctuations. In contrast, 
pairing correlations reduce the number of low-lying lev- 
els (because of the pairing gap) and thus tend to enhance 
the peak height fluctuations. 

There is a regime in the parameter space A/5 — J s /5 



in which the peak spacing distribution exhibits the sig- 
nature of pairing correlations (i.e., it is bimodal), while 
the peak height fluctuations show the signature of ex- 
change correlations (i.e., they are suppressed). This is a 
mesoscopic coexistence regime of pairing and exchange 
correlations. 

(iii) Thermodynamics. Pairing correlations lead to a 
number-parity effect in the heat capacity and in the spin 
susceptibility of the grain. This effect is shifted to lower 
temperatures in the presence of exchange correlations. 
In the fluctuation-dominated regime A/S < 1, the odd- 
even effect is suppressed by exchange correlations. How- 
ever, for A/S > 1, certain signatures of pairing corre- 
lations, such as a peak in the heat capacity for an even 
grain and a local minimum in the spin susceptibility of an 
odd grain, are enhanced by exchange correlations. This 
is explained by the stronger dependence of thermody- 
namic observables on the exchange coupling strength at 
higher temperatures and the much weaker dependence 
at low temperatures. Mesoscopic fluctuations of thermo- 
dynamic observables can wash out the odd-even effects 
for sufficiently small A/S or large J s /S. The fluctuations 
of the spin susceptibility are more sensitive to exchange 
correlations in the fluctuations-dominated regime and be- 
come particularly large as J s /S approaches the Stoncr 
instability. 
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